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;_i I Isospin and U-spin symmetries are useful approximations for studying 

penguin dominance in B meson decays to three kaons, B — > KKK. We 

point out certain subtleties in treating these decays in these approxima- 

^ . tions. Resulting uncertainties are discussed in determining the CP content 

of the final state in i? — > K~^K~Ks, and in relating the CP asymmetry 

I in this process to the value of sin 2/3. 

^ , 

I> . PACS codes: 12.15.Hh, 12.15.Ji, 13.25.Hw, 14.40.Nd 
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Recently the Belle collaboration reported branching ratio measurements for B me- 



■ son decays to three kaon states, K^ K~ K^ , K~^K~Ks, KsKsK^ and KsKsKs P, 

^ . and for B^ decays to K^K^-n^ and K^n^n^ [2]. The BaBar collaboration measured 

\ branching ratios for B^ decays to all three body final states involving charged kaons 

' ' or pions, K^K~K~^, K'^K'ir'^, K^tt^tt~ and n~^7r~^'n'~ JHj. The averaged values of 

the measured branching ratios are, in units of 10~^: 



B{B+^K+K+K~) = 30.8 ±2.1 , B{B^ ^ K+ K^K') = 29.3 ± 5.3 , (1) 

: B{B+ K+KsKs) = 13.4 ± 2.4 , B{B° KsKgKs) = 4.3 ± 1.7 , (2) 

B{B+ ^ K+Ti+Ti-) = 58.4 ±4.4 , ^ 7r+7r+7r~) = 10.9 ± 3.7 , (3) 

B{B+ K+K-TI+) < 6.3 (90% c.l.) . (4) 

These are averaged over the process shown and its CP-conjugate. 

Assuming penguin dominance m B ^ KKK, an isospin analysis was attempted 
by the Belle collaboration [T] in order to isolate the CP-even and CP-odd com- 
ponents of the K^K~Ks final state. This information is useful for studying the 
time-dependent CP asymmetry in this channel [1]. Similar isospin arguments were 
presented subsequently in 0, where a U-spin study relating three body B^ decays 
involving charged kaons and pions was employed in order to estimate deviations from 
penguin dominance. 

In the present Letter we will iterate the isospin analysis for B — > KKK, pointing 
out a subtlety which was overlooked by the above two studies, thereby oversimplifying 
the analysis. It will be shown that these earlier studies made an implicit assumption 
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which goes beyond isospin symmetry. We will argue that, nevertheless, in the pure 
penguin limit an equality holds between the amplitudes of 5° K^K~K^ and 
—>■ K^K^K^ , which is the basis for the CP argument. In order to study deviations 
from penguin dominance in these decays we will employ U-spin considerations which 
were oversimplified in [Hj. We will argue that these deviations, which are partly due 
to electroweak penguin contributions, may be larger than estimated, and introduce 
a sizable uncertainty in the value of sin 2/3 determined from the CP asymmetry in 
K+K-Ks. 

The effective Hamiltonian describing charmless decays B KKK consists of 
operators transforming as a sum of A/ = and A/ = 1. The initial state is pure 
|/ = |), while the final state |/) = \KKK) is a superposition of three isospin states, 
\I{KK) = 0,// = I), \I{KK) = 1, If = 1) and \I{KK) = 1, If = |). Conse- 
quently, these decay processes are described by five independent isospin amplitudes 
[H] corresponding to given KKK momenta. Suppressing the momentum dependence, 

we denote these amplitudes by A^^^f^^'^^ = {I{KK), J/|A/||), and list them as 

0- 1- 0- 1- 1- 

^ , v4q ^ , ^ , ^ and ^ . 

B mesons decay into two kaons and an antikaon in four distinct flavor modes, 
B+ K+K+K-, B° K^K^k", B+ K+K^K° and 5° ^ K+K^K-. The 
first two processes involve identical kaons {K~^ and K^, respectively), while in the 
other two processes the kaons {K^ and K^ in both decays) have different charges. 
The amplitudes for B — > KKK depend on the kaons' momenta. In B^ K^K^K~ 
and in 5+ — > K^K^K one expects different values for amplitudes when exchanging 
the K^ and K^ momenta. Using isospin, one can decompose decay amplitudes for 
B — >■ KKK into isospin amplitudes A^^f^^'^^ describing two kaons and an antikaon 
with given momenta. Thus, one obtains expressions for six decay amplitudes in terms 
of five isospin amplitudes, 

A{K+K+K-),,,,,^ = 2^;'^ - 2aH + a}'^ , (5) 
AiK'K'K^),,,,,, = -2^^'^ - 2aH + aI'^ , (6) 



AiK^K'Tf),,,,,^ = a'/^ - aI''^ - A'i'K AiK a'/^ , (7) 

A{K+K'Tf),,,,,^ = -a'/^ - a'/- + a',''- + a'/^ + a\''- , (8) 

A{K+K'K-),,,,,, = 4U aJ'U A?'^ + + aI'^ , (9) 

AiK+K'K'),^,^,^ = -a'/^ +aI''- -a',''- +a'/^ + a\'^ . (10) 

On the left-hand side amplitudes are specified by the three outgoing particles and by 
their respective momenta. On the right-hand side we have absorbed Clebsch-Gordan 
coefficients in the definition of isospin amplitudes and have suppressed the momentum 
dependence of these amplitudes. 

Let us comment briefiy on Eqs. (j^- lfTUI) . The equal magnitudes and the rela- 
tive signs of contributions of isospin amplitudes in pairs of processes can be easily 
understood from simple considerations. In Eqs. (0) and (jHI) the final states involve 
I{KK) = 1 and are related to each other by an isospin refiection u ^ d, d —u. 
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Consequently, the magnitudes of A/ = and A/ = 1 contributions in these two 
processes are equal and occur with opposite and equal signs, respectively. In Eq. ((7j) 
and (jHl) one interchanges the and momenta in K^K^K . It then fol- 

lows from Bose statistics that contributions to the two amplitudes from I{KK) = 
terms, which are antisymmetric in the isospins of the two kaons, and I{KK) = 1 
terms, which are symmetric in isospin, are equal and have opposite and equal signs, 
respectively. This corresponds to situations in which the orbital angular momentum 
of the KK system in its center of mass frame (which equals the K angular momen- 
tum relative to this center of mass) is odd and even, respectively. The same argument 
applies to Eqs. and pO|) . giving the amplitude for —>■ K^K^K~ in two points 
of phase space where the two kaon momenta are interchanged. 

An interesting consequence of the isospin decomposition is a sum rule between 
B^ and B^ decay amplitudes. The six amplitudes (fSjl- ffTUj) obey one linear relation 
between the sum of three amplitudes for a charged B meson and the sum of three 
amplitudes for a neutral 5, 

A{K+K+K-),,,,,^ + A(ir+i^0Z°),,,,,3 + A{K+ K'^lf = 

AiK'>K'>K'),,,,,, + AiK'^K'K-),,,^,, + AiK-'K'K-)^^^^^ = 3^^ • (11) 

The two sums, in each of which one sums over amplitudes at two points of phase space 
where and momenta are interchanged, are given by the = | amplitude. 
This relation is similar to an isospin relation among the four amplitudes for S"*" and 
S° decays to Kn [7. 

So far our arguments were based purely on isospin symmetry. Let us now study 
the consequences of penguin dominance in 5 — > KKK decays, assuming that the 
dominant term in the AC = 0, AS" = 1 effective Hamiltonian ITU] contributing 
to these decays is a 6 s QCD penguin operator. The assumption of penguin 
dominance implies that one keeps only A/ = terms in Eqs. neglecting 

Ai'^ Ai"^ and A^"^ . Note that this excludes electroweak penguin operators which 
contain a term transforming as A J = 1. We will return to this point when discussing 
deviations from penguin dominance. In the latter approximation one has 

A{K^K^K~),,,,,^ = -A{K'K'Tf ),,,,,, = 2A'/^ , (12) 
AiK-^K'^K'),,,,,,, = -AiK'K-'KX,,,,=Al'^+Ay^ , (13) 

0,- 1,- 

where Aq'^ and Aq^ are antisymmetric and symmetric under interchange of the two 
kaon momenta. The equality of the two pairs of B^ and B^ decay amplitudes follows 
simply from an isospin reflection u ^ d in initial and final states. In this limit the 
two amplitudes involving K^K~ in the final state are, however, different. 

When squaring the amplitudes and integrating over phase space one includes a fac- 
tor I for identical particles in the first pair of processes. The interference between the 

0,- i,i 

two isospin amplitudes Aq'^ and Ag'^ in the second pair of processes, corresponding 
to even and odd angular momenta of the KK system, vanishes. Thus, one finds 

r(S+ ^ K+K+K-) = r(5° K^K^l^) = 2tI'^' , (14) 
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T{B'^ ^ K+K'^K-) = T{B+ ^ K+K°K^) = To' +tI''' , (15) 

0,i l,i 

where Fg'^ and Tq^ are rates corresponding to the two isospin amphtudes. We con- 
clude that, while the two rate equalities p4|) and p5|) follow from penguin dominance, 
not all four rates are equal in this approximation. In particular, the two rates involv- 
ing K^K~ in the final state may be different in general, contrary to arguments made 

in PJIS]. They become equal when Fq'^ = Fq'^ which goes beyond isospin symmetry. 
Experimentally, one has F(5+ ^ K+K+K-)/T{B^ K+K^R-) = 0.98 ± 0.19 

which implies fJ' VrJ' ^ = 1.05^1];^. 

Eq. (fT^ provides the basis for attempting a separation between the CP-even and 
CP-odd components in the final state of the measured process 5° K^K~Ks pp. 
These components correspond to even and odd angular momentum K^K^ states. 
Equal amplitudes in 5° ^ K+K^R- and 5+ ^ K+K^'^ imply that the final 
states in the two processes have the same angular momentum decomposition in terms 
of K^K~ in one process and K^l^ in the other. The probability for a K^T^ being 
in an even angular momentum state, where it decays as KsKs + K^Kl, is given 
by 2T{B+ K+KsKs)/T{B+ K+K^lf). Using Eq. (dS)), this probability is 
given by a ratio of two measured rates 2F(5+ K+KsKs)/T{B^ K+K^K~) = 
1.04 ± 0.20, excluding the (pKs contribution in the denominator 01- This is also the 
probability for a CP-even state in B^ K^K~Ks excluding (pKs- 

The above conclusion, indicating that the final state in B^ K^K^Ks (ex- 
cluding (pKs) is dominantly CP-even, is based on assuming penguin dominance in 
B KKK. The rest of the discussion will address this issue. We will study tests 
for penguin dominance in i? — > KKK, and will evaluate deviations from this approx- 
imation in terms of measurable rates. 

Penguin dominance in 5 — i^vr, suggested in |Hj, was first tested in ^T] by using 
flavor SU(3) and comparing decay rates for B Ktt and B ^ titi. The measured 
rates were also used to estimate the deviation from pure penguin dominance, given 
by a parameter \T'/P'\ ~ 0.2. A similar analysis will be presented here in order to 
relate B KKK to B ^ tttttt and B Ktitt to B ^ KKtt. This will test penguin 
dominance in B ^ KKK. Our arguments differ in detail from those presented in 
m in]. In PP factorization was assumed for three body decays for which no good 
theoretical justification exists, while in j^l a subtlety in using U-spin was overlooked. 
We will also explain the special role of electroweak penguin contributions which were 
ignored in the latter study. 

A useful subgroup of SU(3) permitting relations between B~^ meson decays to 
final states involving charged pions and kaons is U-spin ^2]; under which the pairs 
{d, s), (s, —d), (vr", K~) and {K~^, — vr^) transform like doublets. The AC = 0, AS" = 
1 effective Hamiltonian transforms like a s component {AUs = |) of a U-spin doublet, 
while the AC = 0, AS* = Hamiltonian transforms like a d component = — |) 

of another U-spin doublet. Let us consider the decays of B~^ into the final states 
K^K^K^ , K^TT^K^ , K^T[^Ti~ and tt+tt+tt^. The initial state in these decays is 
pure f/ = 0. The final states from the AU = | transitions are two U-spin doublets, in 
which the two positively charged particles are in f/(++) = and U (++) = 1 states. 
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Therefore, AS = 1 and AS* = amplitudes may be written separately in terms of 
two U-spin amplitudes corresponding to these two states [13]. We will denote AS* = 1 
amplitudes corresponding to f/(++) = and f/(++) = 1 by y4° and A], respectively, 
and analogous AS* = amplitudes by and A\. 

Absorbing Clebsch-Gordan coefficients in the definition of these amplitudes and 
specifying the three outgoing particle momenta, one finds 

AiK^K^K-),,,,,^ = 2Al , (16) 

A(ir+7r+7r-U,,3 = -A^ + A] , (17) 

A{K+n+n-),,,,,^ = A'i + Al , (18) 

A(7r+7r+7r-)p,,,p3 = 2A], , (19) 

AiK+n+K-),,,,,, = A^Al , (20) 

A{K+7r+K-),,,,,, = -A', + Al . (21) 

The similar forms of corresponding AS* = 1 and AS" = amplitudes may be easily 
understood in terms of a simple U-spin refiection d ^ s. The relative signs of terms 
in amplitudes in which two momenta are interchanged follow from Bose statistics. 
Two amplitude relations follow from U-spin, 

A(ir+7r+7r-)p,p,p3 + A(ir+7r+7r-)p,p,p3 = AiK+K+K~)p,p,p, , (22) 
A{K+7i+K')p,p,p, + A{K+n+K')p,p,p, = A{n+n+n-)p,p,p^ . (23) 

On the left-hand-side one sums over amplitudes at two points in phase space where 
the and vr^ momenta are interchanged. 

Squaring Eqs. p6|) - (PT|) and integrating over phase space, one obtains 

r(fi+ ^ K+K+K-) = 2Tl , r(fi+ ^ K+TT+TT-) = r° + , (24) 

r(5+ ^ vr+TT+TT-) = 2ri , T{B+ ^KVK~) = TUT', , (25) 

where r° ^ and ^ are rates corresponding to the two U-spin amplitudes. Contrary 
to arguments presented in 5J, this does not imply an equality between r{B'^ 
K+K+K~) and T{B+ K+tt+tt^), or between T{B+ tt+tt+tt^) and T{B+ 
K^Ti^K~). We conclude that U-spin predictions cannot be tested in simple rate 
equalities. Instead, as shown in |12J, U-spin predicts equal CP rate differences between 
all pairs of U-spin related decays. For instance, the CP rate differences in B^ 
K^K^K~ and B^ tt+tt+tt" are equal in the U-spin symmetry limit. Experimental 
tests of such predictions are quite challenging. 

In order to relate AS" = 1 and AS* = processes to each other, we decompose the 
corresponding effective Hamiltonians into terms multiplying given CKM factors [12|, 

nir = v:,Vuso: + v:,v,sOi , (26) 
kV' = v:,v^,ot+v:,v,,oi , (27) 

where O^'* and O'^'^ are two operators transforming like U-spin doublets. This implies 

A^' = v:,VuA''/ + V:,V,sA'/ , (28) 
AY = + K;K.Ac'' , (29) 
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where A^'^ and ^4°'^ are matrix elements of the U-spin doublet operators Ou and Oc, 
respectively, for U-spin final states with [/(++) = 0, 1. 

Eqs. fOHjl and may be used to test an assumption of penguin and electroweak 
penguin dominance in strangeness changing decays, namely an assumption that the 
second term in Eq. (PHjl dominates the amplitude. We note, however, that in contrast 
to the isospin analysis which assumed QCD penguin dominance, the second term 
includes also contributions from electroweak penguin operators. 

Applying Eqs. ^ and (HH) to S+ K^K+K~ and 5+ tt+tt+tt", respec- 
tively, which involve the amplitudes and Al, and defining momentum-dependent 
amplitudes \P\e'^ = V^VcsAl, \T\e'^ = l/J^Ks^i, one has 

A{B+ K+K-K+) = \T\e'^ + \P\e'^ , (30) 
A{B^ 7r+7r-7r+) = |T|A-^e^^ - |P|Ae*^ , (31) 

where 

^ = = A ^ --^ = 0.226 , (A = 0.22) . (32) 

Vud Vcs 1 - (AV2) 

In the amplitudes for decays to charge conjugate final states the weak phase 7 
appears with opposite signs. Thus, one obtains for the two charge averaged partial 
widths: 

V{B^ ^ K^K^K^) = y"[|T|2 + |P|2 + 2|T||P|cos7cos(5] , (33) 
r(P± ^ TT^TT^TT^) = y[|T|2A-2 + |P|2a2-2|T||P|cos7cos5] , (34) 

where / stands for an integral over three-body phase space. We seek an upper bound 
on 

J \T\V J \P\'\ , (35) 

the tree-to-penguin amplitude ratio averaged over phase space. 

(1) Let us assume that B"^ —>■ tt^ti^tt^ is dominated by the tree amplitude and 
P^ —>■ K^K^K"^ by the penguin. Using the observed value of the ratio 

we then find z = 0.13±0.02 < 0.16 (90% c.L). However, the possibility of tree-penguin 
interference weakens this bound somewhat as we show now. 

(2) The Schwarz inequality implies that in the definition 



J |T||P|cos5 = ^ J |T|2 J |P 



1/2 



(37) 



the magnitude of the parameter ^ cannot exceed 1, and equals 1 when the penguin 
and tree amplitudes are proportional to each other and are relatively real over the 
entire phase space. One can then show that z is a monotonically increasing function 
of ^ cos 7 for any realistic value of TZ. An upper bound ^ cos 7 < 0.74 is based on 
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assuming 7 > 42° ^3]. (The bound on z is not very sensitive to this assumption.) 
Since we do not know b [1^], we shall regard cos 5 as unrestricted. One obtains the 
maximum value of z for ,^cos7 = 0.74: z = 0.19 ± 0.03 or 2; < 0.23 (90% c.L). 

Let us make a few remarks about deviations from penguin dominance which were 
neglected in Eq. (fT5|) : 

1. The bound on the tree-to-penguin amplitude z is based on applying U-spin in 
order to relate K^K^K~ and n^7r^7r~. Large U-spin breaking 
effects are expected to affect a relation between these two processes. This may 
result in a value of z as large as 0.3. 

2. Corrections to penguin and electroweak penguin dominance affecting Eq. (fT^ 
may differ, in both their flavor and momentum dependence, from those esti- 
mated above for 5+ K^K^K~ . Therefore, the bound on z can only be 
used indirectly to set an upper limit of this order on corrections to Eq. (jl3|) 
from tree amplitudes. We conclude that the tree amplitude contribution to 
A{K+K^K-)+A{K^K+^) could be as large as 0.3 of the equal (but opposite 
in sign) penguin amplitudes contributing to these two processes. 

3. Electroweak penguin contributions affect the amplitude equality ()13|) . Such 
terms, which were included in the denominator of Eq. but not in its nu- 
merator, involve the same weak phase as the dominant penguin amplitude. A 
rough estimate of electroweak penguin corrections to the amplitude equality, 
based on Wilson coefficients or on model calculations P, is about 10 — 20%. 

The combined correction to Eq. (|13|) from A/ = 1 tree and electroweak penguin 
amplitudes is hard to calculate, and depends on the interference between tree and 
electroweak penguin amplitudes which may be constructive except in special cases 
[T^ . Constructive interference could imply an overall correction of 40 — 50% in the 
most pessimistic case. 

The effects of electroweak penguin and tree amplitudes on the time dependent 
CP asymmetry in B^ —>■ K^K^Ks are of two kinds. Both amplitudes affect the CP 
structure of the final state / = K^K~Ks, thereby multiplying by a dilution factor 
the coefficient Sf of the sin Ami term, where Sf = sin 2/3 for a purely CP-even state 
excluding (pKs- One way of determining the dilution factor is through a partial wave 
analysis in the angular momentum of the K^K~ system. This may separate even 
and odd angular momenta corresponding to even and odd CP states. 

The tree amplitude (Tj) in B^ — > K" ^ which has a weak phase different 
from that of the penguin and electroweak penguin amplitudes (combined to P/), has 
another effect. [We use the same convention for Tf and Pf as in (jHUj) ]. It modifies 
5*/ from a value sin 2/? for a pure CP-even state (excluding (pKs) to sin2/5efT, (3^^ = 
(3 + 5(3, and introduces in the asymmetry a cos Ami term with coefficient Cf. These 
corrections depend on the tree-to-penguin ratio of amplitudes and on the relative 
strong phase between these amplitudes (5/), both varying over phase space. We 
define 




(38) 
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Keeping only linear terms in Zf, one obtains for Cf and 6 [3 expressions which 
generalize those obtained for two body decays [T7j 

Cf ~ 2zf^ssm^ , (41) 
6p ~ Zf^csin-f . (42) 

We estimated that Zf may be as large as 0.3. The largest possible effect on Sf occurs 
when l^cl = 1; corresponding to penguin and tree amplitudes which are proportional 
to each other and are relaively real over the entire phase space. In this case \6(3\ could 
be as large as 17°, and Sf could lie in the range 0.2 — 1. The Schwarz inequality, 
+ — 1) implies = when |^c| = 1, hence Cf = 0. The other extreme 
but unlikely case, of a maximal |C/| = 2zjsin7 and a minimal 6f3 = 0, occurs 
when C,s = ±1, = 0. This corresponds to penguin and tree amplitudes which are 
proportional to each other and are out of phase by 7r/2 over the entire phase space. 

M. G. wishes to thank the CERN Theory Division and the Enrico Fermi Institute 
at the University of Chicago for their kind hospitality. This work was supported 
in part by the United States Department of Energy through Grant No. DE FG02 
90ER40560. 
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